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ABSTRACT 



Fjortoft's graphical method of prognosis of the 500-mb 
surface assumed relative vorticity change as being due to 
advection only. This study considers the total 24-hour -rela- 
tive vorticity change as consisting of two components, that 
due to 100# advection and that due to non-advective changes. 
Objective methods are developed for the prognosis of the 
latter. These changes are then added to the advective changes 
to give the total 24-hour relative vorticity change from which 
the 24-hour height change can be recovered. 

The 500-mb velocity field is subjectively categorized 
into five types. For each type the non-advective relative 
vorticity change is correlated with a parameter determined 
from the spreading of the contours. Results for the five 
categories, when applied to data for the months of February 
and March 1957, are as follows: 







Correlation 


Case 


Velocity Field 


Coefficient 


I 


Marked Cyclonic Curvature 


0.60 


II 


Moderate Cyclonic Curvature 


0.71 

/ 


III 


Negligible Curvature 


0.79 


IV 


Moderate Anticyclonic Curvature 


0.56 


V 


Marked Anticyclonic Curvature 


0.10 



Quantitative objective forecast rules are developed from the 
regression line’s relating the non-advective relative vorticity 
change and the parameters for each of the five cases. 

Using the developed technique the results of five 24-hour 
500-mb test forecasts for the United States are presented. 
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These tests show an Improvement of four feet In forecast 
height value for the entire forecast area as compared to 
the Fjortoft method. However, in considering only those 
points at which the computation was made an average improve- 
ment of 34 feet is noted. 

Finally the persistence of non-advective relative 
vorticity changes is considered for its value as a height- 
change forecast tool. The results of seven forecasts pre- 
pared using this assumption indicate an improvement of 32 
feet over the Fjortoft method for the entire forecast area. 

The writer wishes to express appreciation to Professor 
W. D. Duthie for guidanoe and assistance in the preparation 
of this paper. 
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CHAPTER I 



INTRODUCTION 

Fjortoft^lJ developed a graphical method of fore- 
casting the 500-mb constant pressure surface by the follow- 
ing procedure: 1) determine the distribution of relative 

vorticity at the 500-mb level'; 2) advect this vorticity with 
the mean 500-mb flow; 3) obtain the 24-hour vorticity change 
by subtraction of these two fields; 4) recover the 24-hour 
height-change field from the 24-hour vorticity-ohange field; 
and 5) arrive at the prognostic chart by adding the height- 
change field to the current 500-mb chart. The assumptions 
involved in this procedure include barotropy, the geostrophic 
approximation, the neglect of the change of coriolis para- 
meter and map scale with change in latitude, and the conser- 
vation of the mean 500-mb flow during the forecast period. 

Many less restrictive models have been developed. 

Fjortoft, in the paper cited above, compensated for the change 
of coriolis force and map scale with latitude by adding a 
field, which is a 'function of latitude only, to the relative 
vorticity field and advecting this quantity. In another paper, 
Fjortoft j^2^has developed a baroclinic model that assumes a 
functional relationship in the wind profile. Hesse [^3^j has 
included baroclinity with vertical motion due to terrain. 
Without question these more sophisticated graphical tech- 
niques will prove superior to arbitrarily seleoted statisti- 
cal parameters. However, for the meteoroligist in the field. 
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the time involved in the preparation of the routine fore- 
cast must be held within strict limits. It is felt that 
techniques more complicated than the simple Fjortoft method 
would not be within the capabilities of the average Naval 
Weather Station. 

The Air Weather Service jjtj has tested the simple 
Fjortoft method as a means of objectively forecasting the 
500-mb contour field for a 24-hour period. In the test 
consisting of 541 samples the Fjortoft method gave an 
average absolute error of 208 feet in the forecast 500-mb 
height, whereas subjective techniques employed at the TJSAF 
Weather Central gave an average error of 207 feet. It was 
concluded that: 

objective use of the graphical method presented 
.by R. Fjortoft produces 24-hour 500-mb prognostic 
charts of about the same quality as those produced 
by skilled forecasters in a weather central. 

It was further concluded that the time required for the prepa- 
ration of the complete forecast was within the capabilities 

/ • 

of the average weather detachment. 

Because of its assumptions of barotropy, the Fjortoft 
method is incapable of forecasting the changes in the in- 
tensities of circulation; henoe it will not show the rapid 
and spectacular deepening of lows aloft, which are ascribed 
to the baroclinic state of the atmosphere. This paper is an 
attempt to develop statistical methods which will indioate 
the magnitude of baroollnlc changes in the upper air and, 
together with the Fjortoft procedures, provide an improved 
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objective forecast routine for the meteorologist in the 
field. 
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CHAPTER II 



THEORETICAL DEVELOPMENT 



The vorticity equation is expressed in isobaric 
coordinates as follows: 



dr 



~(J \f)v-v 




where / is the vertical component of ^he relative vorticity, 
f is the coriolis parameter, V is the horizontal velocity, 
p is the pressure, and co - ^ . Vectors are indicated by 
overlining, and the vector operator V 1 u + Sy 1 ^ ^ 

has its usual vector meaning. Expanding the left hand side 
of this equation. 



*M£lQ + j/' 7 1 j i 

dt dr 

= it + v- %J + v- * CO 

dt . H H dp-' 




and combining with the right hand side we have: 

- - V- V„ i ~ v f -60 - tti /) 7 ' / - VcO * • k (3) 

or the looal rate of change of relative vorticity is equal 
to the sum of terms representing horizontal adveotion of 
relative vorticity, horizontal adveotion of the change of 
ooriolis parameter with latitude, vertical adveotion of 
absolute vorticity, generation of absolute vorticity due to 
divergence in the velocity field, and generation of absolute 
vorticity by the "tilting' 1 term. If we let the sum of these 
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last four terms be equal to R, this equation may be ex- 
pressed as 



FJortoft, In his simplest model, considered R negligible 
and forecast the local rate of ohange of relative vortloity 
by advecting the relative vorticity at the mean geostrophlc 
wind velooity. If statistical methods were utilized to 
determine the quantity R and this value added to the advected 
field, the forecast would be improved. 

The purpose then is to develop an objective prognostic 
routine consisting of the following steps: 1) determine the 

relative vorticity field; 2) advect this field with the mean 
geostrophlc wind; 3) graphically calculate the change in 
relative vorticity due to adveotion; 4) add to this change 
field an empirically determined correction to compensate for 
' the neglected terms in the vorticity equation; and 6) recover 
the height-change field and the contour field from the 
corrected vorticity-ohange field. 

Well known methods {[ 5 J are available for graphically 
computing the field of relative vorticity. Making use of 
the geostrophlc approximation, the relative vorticity can 
be expressed as 







i-7 U 

where Z is the contour height and \J is the Laplaoian 
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operator in the plane. Written in finite difference form 

^ ' id'- ^ (&) 

we have a means for determining the field graphically, where 
the subscripts denote values of Z at the four sides and at 
the center of a square grid, and d is the grid distance. If 
the small variation of f with latitude is Ignored, the rela- 

V 

tive vorticity at any point may be represented by 

J ■- 2 -■? ( 7 ) 



where 






is the space-mean contour height. 

Now combining equations (4) and (7) 
by the more conservative velocity field 







replacing 
we have 



V 




where is the geos trophic mean wind. 

If the equation for relative vorticity (5) be differen- 
tiated with respect to time, we have 




The solution to this partial differential equation offers 
the means of recovering the height-change field from the 
forecast vorticity-change field. Fjortoft £V] has developed 
a graphical solution to this equation by expanding in a con- 
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verging infinite series and then neglecting the smaller 
terms. If all but the first two terms are neglected, the 
height change is approximated by 

6? = - v/ -a&t 

where ^ denotes the 24-hour change and the bar denotes the 
space-averaged value of the quantity. For a less accurate 
approximation, all but the first term may be ignored giving 

A 2 = -A/ (JJ) 

or the 24-hour height change is equal to the negative of the 
24-hour vorticity change. 



t 
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CHAPTER III 



PROCEDURE 

Used in this research was a series of 500-mb charts of 
the United States and the Pacific Ocean prepared by the staff 
analysts at the U. S. Naval Postgraduate School, Monterey, 
California, during the months of February and March 1957. 

The field of 24-hour vorticity change due to advection 
was determined by the Fjortoft method. This field was then 

s'* 

subtracted from the actual vorticity-change field as calcu- 
lated from the next day*s actual vorticity. This field, R 
in equation (8), was considered to be the vorticity change 
due to baroclinic influences in the atmosphere. The fore- 
casting of this quantity will be considered in the next 
chapter. 

In constructing vorticity charts by graphical means 
a grid distance of 6 degrees of latitude at 60 degrees North 
was considered to give best results. Isopleths of relative 
vorticity were advected at 100 percent of the space-mean 

T 

500-mb wind. Most authors {^4, 6 ] suggest 80 percent as 
being best in order to compensate for the use of the 500-mb 
level as the level of non-divergence. However, as the objective 
of this paper is to forecast the effects of divergence on the 
vorticity field, 100 percent of the mean wind was used. 

Some difficulty was encountered in the graphical tech- 
niques of adding and subtracting due to fields consisting of 
many small, closed centers. This difficulty was usually 
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overcome by drawing intermediate iaopleths. However, on a 
few occasions the field was still not adequately determined; 
but it was evident that this occurred only when the vorticity 
change was small and due entirely to advection. 
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CHAPTER IV 



PORECASTINO VORTICITY CHANGES 

Aa discussed in Chapter II, two of the sources of non- 
adveotlve relative vorticity changes are divergence in the 
velooity field and vertical motions in the atmosphere. An 
attempt will be made to forecast these changes by purely 
statistical means. 

Scherhag £ 7 J states that where upper contours converge, 
an "entranoe zone”, pressure rises will occur due to trans- 
verse convergence, and where upper contours diverge, a "delta 
zone", pressure falls will occur due to transverse divergence. 
CMConnor [^8] states that convergence and upper height rises 
are associated with low speed winds approaching straight or 
cyclonically curved strong contour gradients or with high 
speed winds approaching anticyc Ionic ally curved weak contour 
gradients, and that divergence and upper height falls are 
associated with high speed winds approaching cyclonically 
curved weak contour gradients. These two statements provide 
the parameters that will be correlated with non-adveotive 
vorticity changes in order to forecast the magnitude of the 
changes due to divergence in the velocity field. 

0*Neill "} conducted an investigation as to the validity 
of the above statements for forecasting height-change centers 
and concluded that, although qualitatively correct, the rules 
devised did not show sufficient correlation between observed 
and forecast height change to provide a quantitative fore- 
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cast. However, It is felt that even though no high corre- 
lation existed between total height-ohange and selected 
parameters, such a correlation may exist when non-advective 
vorticity change is the dependent variable. 

An attempt was made to measure the radius of curvature 
as one of the parameters to be used. Because of difficulties 
in determining where to make the measurement this was abandon- 
ed in favor of using five subjeotlve classifications: marked 

cyclonic curvature, moderate cyclonic curvature, negligible 
curvature or straight contours, moderate anticyclonic curva- 
ture, and marked anticyclonic curvature. These five cases 
will be examined in detail. 

Case I - Marked cyclonic curvature . There were a total 
of 21 configurations that were considered to fall within 
this classification in the series of map3 examined. All 
were associated with major troughs, one located over the 
west coast of the United States moving slowly eastward and 
the other over the Western Atlantic which moved off the chart. 

The parameter used to measure the divergence of contours 
is defined as xi s d^ - d^ where dt is the distance in 
degrees latitude between 400 feet of contour at the trough 
line, and d^ is the distance between 400 feet of contour at 
the upstream inflection point of the same contours. The 
contours selected are those in the geometric center of the 
flow around the trough. Figure 1 shows the configurations 
involved and the method of measuring x^. 
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Multiple correlation techniques of R with two independent 
variables, magnitude of the vorticity change and location 
relative to the trough line of the vorticity -change center, 
were attempted, but no significant correlation was found to 
exist^* j n 19 cases there were closed centers of non-adveciJve 
vorticity change, the mean position of which was located 
approximately on the middle one of the three involved contours 
at a distance downstream from the trough line of five degrees 
of latitude. Therefore, this was taken as the position at 
which R was measured. 

Figure 2 is the scatter diagram of the values obtained 
in the investigation of Case I examples with the regression 
line of R on x^, R s 116 -f 46 xj_. The correlation coefficient 
between the two variables was calculated to be 0.60. 

It is to be noted that there were only 4 cases in which 
the value of x^ was negative, that is, the geos trophic wind 
at the trough was greater than at the inflection point. If 
the four cases were eliminated, the correlation coefficient 
would be improved. 

Case II - Moderate cyclonic curvature . Fifteen examples 
were observed in this classification. These were areas of 
large longitudinal extent with generally smooth flow and 
little change in contour spacing* 

^These techniques were applied in all cases with similar 
results. 
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The values of the parameter selected for the measurement 
of divergence, x 2 , were obtained as in Case I. There were 
only four closed centers of non-advective vorticity change 
associated with this configuration of contours, the mean 
position of which was different from that in Case I; however, 
in the interest of consistency, the 3ame location for measure- 
ing R was used in both cyclonic cases. 

Figure 3 is the scatter diagram for the 15 Case II ex- 
amples, A regression line of R on x 2 , R s 24 -f 78 Xg, was 
obtained. These data gave a correlation coefficient of 0,71. 

Case III - Negligible curvature . There were 18 examples 
observed under this classification. They were located in 
broad troughs or ridges and in areas of extensive zonal flow. 

The value of the parameter X 3 is defined as X 5 s dg - d-^ 
where d^ is the distance in degrees latitude between 400 feet 
of contour at the upstream end of the straight contours, and 
dg is the distance in degrees latitude between the same con- 
tours at the downstream end. The contours used were those in 
the center of the flow. The value of R was measured at the/ 
downstream end of the center contour. Figure 4 indicates the 
parameters used. 

The scatter diagram, figure 5, gives the results of the 
Case III examples with the regression line of R on X 5 , 

R - 19 + 45 X 3 , with a correlation coefficient of 0.79. 

Case IV - Moderate anticyclonic curvature . Ther^ were 

15 examples in this category, located in areas of broad zonal 
flow. There was little change in contour spacing involved 
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in these areas 



The value of the divergence parameter, X4, is defined 
as m dp - d^, where dp is the distance in degrees latitude 
between 400 feet of contour at the ridge line, and d^ is the 
distance in degrees latitude between the same 400 feet of 
contour at the inflection point upstream. The value of R 
was taken five degrees of latitude downstream from the ridge 
line on the contour of least height value of the three used. 
Figure 6 indicates the method of measuring the parameters. 

The scatter diagram, figure 7, indicates the results of 
the Case IV investigation with the regression line' of 
R on X4, R ■ -40 + 53x4* A correlation coefficient of 0.56 
was obtained. 

Case V - Marked antlcyclonlc curvature . There were 24 
configurations that were classified in this category, all 
directly upstream or downstream from the major troughs of 
Case I. 

The value of the parameter X5 was obtained in the same 
manner as X4 in Case IV. There were 17 closed centers of 
non-advective vorticity change involved with a mean position 
5 degrees of latitude downstream from the ridge line near the 
contour of least height value. R was measured at this point. 

Figure 8 Is the scatter diagram for Case V with the re- 
gression line of R on X5, R ■ -155 4- 6x5. A correlation 
coefficient of 0.10 was obtained. 

The results of the above part of the investigation are 
summarized in Table I. 
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Case I 



Number of Examples 
Parameter (°Lat.) 

Regression Line 
Correlation Coefficient 
Mean value of R (feet) 

Standard deviation of R (feet) 

Case II 

Number of Examples 
Parameter (°Lat.) 

Regression Line 

Correlation Coefficient 
Mean value of R (feet) 

Standard deviation of R (feet) 

Case III 

Number of Examples 
Parameter (°Lat.) 

Regression Line 
Correlation Coefficient 
Mean value of R (feet) 

Standard deviation of R (feet) 

Case IV 



Number of Examples 
Parameter (°Lat.) 

Regression Line 
Correlation Coefficient 
Mean value of R (feet) 

Standard deviation of R (feet) 

Case V 



Number of Examples 
Parameter (°Lat.) 

Regression Line 
Correlation Coefficient 
Mean value of R (feet) 

Standard deviation of R (feet) 



21 

*L • d t - di 
R = 116 -+- 46x x 

.60 

165 

153 



15 

x 2 - dt - di 

R a 24 + 78x 2 
.71 
20 
95 



18 

X3 - d 2 - di 
R a 19+45x3 
.79 
30 
50 



is 

X4 a d r - d i 

R a “40 ~h 53X4 
.56 
-30 
80 



24 

X5 m d r - di 
R a -155 + 6x5 
.10 
-150 
130 



TABLE I 
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As Indicated earlier, in addition to the divergence in 
the velocity field, there is another souroe of change in the 
local relative vorticity that would appear to lend itself to 
statistical methods: that due to vertical motions in the 

atmosphere* An attempt was made to derive a correction for 
the effects of vertical motions* 

Case I samples were divided into two categories: Those 

in which the major trough was on the lee side of the conti- 
nental divide and those in which the major trough was on the 
windward side* Lines of regression and correlation coefficients 
of R were computed for each of these two categories* There 
were no significant differences* However, the number in each 
case was very small* Results would be more conclusive if 
examples were taken of the two categories in Case III, the 
straight flow case, due to its higher correlation factor, 
but the series of maps used in the research was not adaptable 
to such analysis. Time considerations prohibited further 
testing. 

Conclusions * The correlation coefficients in Cases I, 

II, and III are significant at the 95$ level* It would be 
indicated that in the regions of cyclonic ally curved or 
straight contours a corrective field can be added to the 
advected field of relative vorticity that will improve the 
forecast for a 24-hour period. The correlation coefficient 
in Case IV Is significant at the 90$ level, and in Case V 
there was no statistical significance* In the cases of anti- 
cyolonically curved contours, the radius of curvature of the 
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contour appears to be much more critical, but with moderate 
curvature a correction oan be applied that will improve the 
forecast. With marked curvature a parameter that more accu- 
rately determines the deviation from the gradient wind will 
have to be used. 

It should be noted that the local rate of change of 
relative vorticity in the absence of marked curvature of 
contours is due primarily to advection. Table I shows the 
small values of the means and standard deviations of non- 
advective changes in these cases. 



25 



I 

I 



I 



CHAPTER V 



INDEPENDENT TEST RESULTS 

The 500-mb prognostic techniques developed in Chapter IV 
were tested with a series of five independent forecasts pre- 
pared during the month of April 1957. 

The 24-hour relative vorticity change due to advection 
was determined by the Pjortoft method using 100# of the mean 
500-mb wind as the advective field. The 24-hour relative 
vorticity change due to other sources was prepared by plotting 
the values of R determined from the regression formulas, and 
then subjectively drawing the R field over the entire map. 

i 

_ S 

The sum of these two fields was considered the total 24-hour 
relative vortioity-change field, from which, the 24-hour 
height-change field was recovered by the application of equa- 
tion (10). 

The prognostic charts prepared were verified on the basis 
of forecast and observed 24-hour height changes. Thirty-six 
grid points on each map, for a total of 180 cases, were used 
in this phase. As a basis of comparison of the value of the 
method of prognosis, forecasts were also prepared by the un- 
modified Pjortoft method. 

The root-mean-square errors for the two methods showed 
no significant difference when all the 180 points were con- 
sidered, being 146 feet for the Pjortoft method and 142 feet 
for the other. However, if only those points were considered 
for which a value of R had been actually calculated, a signifi- 
cant improvement was made. There were a total of 37 such 
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points for the five maps. The root-mean-square error in 
forecast height-change for these points was 135 feet as com- 
pared to an error of 169 feet for the forecast prepared by 
the unmodified FJ or toft method. 

Conclusions . The results of the test forecasts indicate 
that the objective prognosis of the 500-mb const suit pressure 
chart cannot be accomplished by the methods developed herein 
for the following reasons: 1) subjectivity enters into the 

classification of the contour pattern; 2) the rules are appli- 
cable at only a few points; 3) it is impossible to draw an 
extended area chart of the values of R given the values at 
only these few points; and 4) the method fails in areas of 
marked anticyclonic flow. However, for limited application, 
it is believed that the values determined from the configu- 
ration of the contour field can be subjectively employed to 
Improve the Fjortoft prognostic chart. 
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CHAPTER VI 



FURTHER STUDIES 

A further study was conducted in order to attempt to 
discover continuously variable parameters for forecasting 
the non-advec tive relative vorticity change* Used in this 
phase of the research was a series of maps for November 1955* 

There was no indication of any continuity in the paths 
of the centers of values of R. Hence extrapolation was of 
no prognostic value. 

A measure of the baroclinity of the atmosphere is the 

number of intersections of isotherms and contours on a con- 
« 

stant pressure surface. An attempt was made to use this 
principle for forecasting in two ways. First, the actual 
density of the solenoids was plotted against the values of R. 
This was done in several different ways, for example, by 
counting the number of intersections of isotherms and contours 
within a radius of three degrees of latitude of the point at 
which R was measured. Second, multiple correlation was 
attempted between the gradients of temperature and height 
value and R. Neither of these attempts met with success. 

Persistence of the values of R was determined to be the 
most successful prognostic tool in the forecasting of R. 

That is to say, for forecasting the value of R at any point 
on the chart, best results were obtained by making the 



assumption that 



JR . 



Jt 
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For the seven maps considered, using 100# of the mean 
flow for the advection of relative vorticity, the Fjortoft 
method gave a root-mean- square error of 211 feet in forecast 
height change as approximated by equation (11). Using per- 
sistence as the forecast of R this error was reduced to 179 
feet. 

As mentioned in Chapter I, Fjortoft £l] indicated a method 
of compensating for the change of coriolis force and map scale 
with latitude. In this improved model, equation (8) was re- 
placed by 

^.(2-2 +J( /P)] * -(Vq +£ 

where a v . 2 . ■ . . / 

<W) -f 

rj y 0 xrg 

and is the earth's rotational velocity, d is the grid 
distance, m is the map scale constant, Cf is the latitude, 
g is gravity, and 

^ VJ (?) 

is a steady easterly wind. Also mentioned earlier was the 
evidence that 80# of the advective field gives better results. 

The same forecasts were prepared using this more refined 
Fjortoft method and the recommended 80# wind. The root-mean- 
square error in forecast height change was 203 feet. Using 
persistence as the forecast for R, this was reduced to 177 
feet. 

It appeared that persistence of R relative to the contour 
field, rather than to the actual geographic position, might 
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give superior results, but evidence was inconclusive. 

Conclusions , Any conclusions reached concerning the 
above results must be qualified ones because of the meager- 
ness of the tests on which they are based. Although the 
improvement in root-mean-square error in the forecast height 
change is statistically slgnif leant , only one series of maps 
was studied, A different synoptic situation may conceivably 
produce entirely different results. Further testing will have 
to be conducted before objective application of the rules 
would be warranted. 

The primary effort in this research was directed at 
developing forecast rules based on the configuration of the 
contour field. The methods developed failed in objective 
application for the reasons given in Chapter V, As a result 
of the further study, these objections have been overoome by 
assuming the persistency of R, Preliminary results indicate 
that this assumption is a valid one and can be objectively 
applied to prognosis of the 500-mb contours. Final con- 
clusions can be reaohed only after further testing. 
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